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Abstract 

The superconducting phase of the 2D one-band Hubbard model is studied 

within the FLEX approximation and by using an EHashberg theory. We 

investigate the doping dependence of Tc, of the gap function A(k, lv) and 

of the effective pairing interaction. Thus we find that becomes maximal 

for 13 % doping. In overdoped systems decreases due to the weakening of 

the antiferromagnetic correlations, while in the underdoped systems due to 

the decreasing quasi particle lifetimes. Furthermore, we find shadow states 

below Tc which affect the electronic excitation spectrum and lead to fine 

structure in photoemission experiments. 
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There is still no definite understanding of the pairing interaction in the High-Tc super- 
conductors. However, improved experimental techniques like phase sensitive measurements 
of the superconducting order parametei0 and angular resolved photoemission (ARPES) i 
have revealed important new structures in the electronic excitation spectrum that might 
help to clarify the nature of superconductivity in the cuprates. The observation of shad- 
ows of the Fermi surface (FS) in the normal state by Aebi et a/.i and the strong evidence 
for a dnc2_y2 pairing symmetryEUl indicate that there might be a strong interdependence 
between the occurrence of high transition temperatures and antiferromagnetic spin fluctu- 
ations. Furthermore, the interesting dip structures in ARPES spectra of Bi2Sr2CaCu208+(5 
(Bi2212) found for the first time by Dessau et al% and confirmed by other measurementsili'i 
is currently debated and could be a further fingerprint of the pairing mechanism. 

From a theoretical point of view the two dimensional one-band Hubbard Hamiltonian 
serves as a basic model to describe the highly correlated electrons within the Cu02 planes. 
The first important steps to investigate the superconducting properties of the cuprates within 
this approach was done simultaneously by Monthoux et a/.0 and Pao et a/.i, who found a 
dx2_y2 pairing state for low temperatures within a pertubative ansatz, namely the fluctua- 
tion exchange approximation (FLEX). Despite these interesting results the superconducting 
excitation spectrum and in particular its doping dependence are far from being understood. 
Moreover, it is still unclear if this microscopic approach can explain the observed optimal 
doping concentration in the High-T^ superconductors. 

In this article we present new results for the doping dependence of the superconducting 
state of the Hubbard Hamiltonian. We discuss the resulting phase diagram and find that 
Tc becomes maximal for a doping concentration of 13 %. By determining the excitation 
spectrum directly on the real frequency axis we demonstrate that shadow bands as proposed 
by Kampf and Schriefferi for the normal phase are also present below Tc. We show that 
they have an important impact on the pairing mechanism, on Tc and on the gap function 
A(k, In addition, they cause observable effects in the spectral density of states and are 
responsible for fine structures in ARPES experiments. 



Our theory is based on a strong coupling Eliashberg approach for the one-band Hub- 
bard Hamiltonian with nearest neighbor hopping integral t = 0.25 eV, bare dispersion 
eo(k) = {—2t[cos{kx) + cos(A;j;)] — /i) with chemical potential fi and local Coulomb repulsion 
U = At. The superconducting state is treated in the Nambu formalism, where the Greens 
function G(k, cu) is a 2 x 2 matrix that can be expanded in terms of Pauli matrices. The cor- 
responding expansion coefficients of the diagonal electronic self energy are — Z(k, u)) and 
x(k, cu), whereas 0(k, cj) = A(k, C(j)Z(k, cu) is the coefficient of the anomalous off-diagonal 
self energy that vanishes above Tc and A(k, cj) is the gap function. These three functions 
were calculated self consistently by using a newly developed numerical method for treating 
the FLEX equations directly on the real frequency axis. The calculations were performed 
on a (64 x 64) square lattice in momentum space and by using 4096 equally spaced energy 
points in the interval [—30 t, 30 t] leading to a low energy resolution of 0.014 t ^ 4 me 

In the Fig. 1(a) we present the doping dependence of Tc, where x is the doping con- 
centration and n = 1 — X the occupation number per site. Note that the order parameter 
was found to have a da.2_y2 symmetry for all doping values. For large x, increases up 
to Tc = 79 K aX X = 0.13, whereas below the optimal doping starts to decrease for 
X < 0.13. Consequently, three different doping regimes with qualitatively different behavior 
occur within the Hubbard model: The underdoped {x < 0.13), the optimally doped and the 
overdoped {x > 0.13) compounds. In the following we will demonstrate (i) that the optimally 
doped systems are characterized by a constructive interference of superconducting and an- 
tiferromagnetic excitations, resulting from the comparable time scales of both phenomena, 
(ii) that the suppression of Tc in the overdoped systems is due to the increasing weakening 
of the effective pairing interaction, and {Hi) that in the underdoped systems the dynamical 
character of the spin fluctuations becomes less pronounced causing Tc to decrease. Here the 
rather stiff arrangement of antiferromagnetically correlated spins and their interplay with 
the (i-wave pairing state enhances the formation of shadow states, i.e. states resulting from 
the coupling of wave vectors k at the FS and its shadow at k + Q (Q = (vr, 7r))i. 

Important information about the origin of the doping dependence of Tc and the optimal 
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doping can be obtained from the effective interaction Vs(k,u!) below TcQ. In Fig. 1 (b) 
we compare Kj(k, tu) with y„(k, cj) that was obtained at the same temperature, but with 
= 0, for cij = and k points near Q where K(k, 0) is maximal. For non-optimal doping 
we find that the interaction and in turn the antiferromagnetic correlations increase only 
slightly below T^, whereas for compounds close to x = 0.13 an enhancement up to 40 % 
occurs with respect to the normal state. It is interesting that this happens although our 
calculation yields that K(k, cu) increases continuously for decreasing doping even below the 
optimal doping. The enhancement of the pairing interaction is consistent with the interesting 
doping dependence of Aq = A(T = 0) = 0(T = 0)/Z{T = 0)0. The superconducting gap 
has a maximum slightly above the optimal doping and decreases sharply in the underdoped 
compounds, because the quasi particle scattering increases for smaller doping due to the 
strong antiferromagnetic correlations via Z, whereas does not change significantly for 
X < 0.13. In addition the 2Ao curve intersects the curve at x = 0.13, where r is 
the lifetime of the quasi particle, and r^^ = Im S(k, 0) with FS-momentum and 
self energy S(k, 0) = uj{l — Z(k, cu)) + x(k, cu) at T^. In Fig. 1 (b) one sees that 
increases monotonously with decreasing doping and saturates for small and large values of x 
as observed in transport measurements^. Therefore the occurrence of an optimal doping in 
our results is related to an interference of the typical lifetime of the Cooper pairs ~ Ag ^ and 
the lifetime r of the quasi particles: First increases for decreasing doping, since Vs{]<.,uj) 
increases and since r is large enough to guarantee an effective Cooper pair formation. Then 
below the optimal doping, decreases again even so K(k)^^) continues to increase, since r 
becomes too small and there are no well defined quasi particles during the pairing process. 
Notice, that there are also remarkable variations of the ratio 2Ao //csTc upon doping shown 
in Fig. 1(c). The deviations of 2Aq/ ksTc from the BCS value demonstrate the qualitative 
new changes of the superconducting state for decreasing temperatures, which are mostly 
pronounced at the optimal doping, where a strong increase of the pairing interaction occurs. 

The interdependence of dynamical antiferromagnetism and superconductivity leads not 
only to an optimal doping, but also to a different behavior of the gap function A(k, cj) in 
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overdoped and underdoped compounds. In Fig. 2 we plot A(k, u) for different frequencies 
along a path in the Brillouin zone as indicated in the insets. By comparing our results 
with the ansatz Ad ~ cos{kx) — cos{ky) one sees deviations from the simple d-wave behavior. 
These are mostly pronounced for small frequencies and in particular for underdoped systems, 
where the gap function is enhanced and largest at the FS and at its shadow. By transforming 
our results for A(k, c<j) to real coordinate space, we find that this corresponds to a pairing 
dominated by nearest-neighbor interactions in the case of x = 0.16. However, the appear- 
ance of higher harmonics in A(k, cj) for x = 0.09 indicates that pairing processes between 
more distant antiferromagnetically correlated spins come into play. This results from the 
increasing stiffness of the nearest-neighbor spin correlations, which do not contribute any 
longer as strongly to the spin fluctuation mediated pairing interaction. 

Further insight into the effect of the short range antiferromagnetic order on the excita- 
tions can be obtained by investigating the frequency dependence of the quasi particle decay. 
This is dominated by the coupling of states at the FS and its shadow. By inverting the ma- 
trix Greens function G(k, u) one can define an effective electronic self energy that includes 
the off-diagonal contributions and is equal to S(k, cu) in the normal state0: 

E^(k, u) = coil - Z(k, u)) + x(k, u;) + ^ Jt^^u^^^^^ ^ (I. ^ " 

In Fig. 3 we compare the imaginary part of E^(k, u;) at the FS and at its shadow for an 
overdoped and an underdoped compound above and below T^. In the normal state we find 
for X = 0.16 that Im S(^(k, cj) is linear in u down to energies ~ 8 meV ~ 80 K, which 
is consistent with t he marginal Fermi-liquid (MFL) theorylli. For x = 0.09 we find besides 
an overall increase of the scattering rates a minimum (III) at the shadow at the FS that is 
a precursor of the singular behavior of ^^(k, cu) in the antiferromagnetic state. Thus the 
underdoped compounds are clearly not in agreement with the properties of an MFL. Below 
Tc and for x = 0.16 a sharp minimum (I) becomes visible at the Fermi energy {uj = 0) 
when k is close to the FS which is the analogue of the corresponding 5-function in the 
BCS theoryi. Thus minimum I shifts well below the Fermi energy when k approaches the 



shadow of the FS. Furthermore, its intensity decreases and becomes almost equal to the 
additional minimum (II) at u ^ +50 meV in Fig 3(a) and (b). Note that minimum II 
vanishes for k points away from the FS and its shadow (not shown) and is purely caused 
by the opening of the superconducting gap and a related shift of spectral weight in the 
effective interaction. Hence, this dip reflects the enhancement of the antiferromagnetic 
coupling below Tc as already discussed in the context of Fig. 1(b). A similar feedback effect 
between spin fluctuations and d-wave pairing state occurs in the underdoped system. Here, 
the strongly enhanced scattering rates and the smaller superconducting gap lead to a rather 
small minimum I at the FS with respect to a; = 0.16 and to a much weaker minimum II above 
the Fermi energy, whose energy shifts to the Fermi level for decreasing doping. However and 
more interestingly, the remarkable enhancement of minimum III shows that it is not a simple 
superposition of the precursor of the antiferromagnetic singularity and the superconducting 
singularity, but caused by a real interdependence between both phenomena. 

To demonstrate that the interplay between dynamical antiferromagnetism that causes 
shadow states in the cuprates and superconductivity has not only consequences for the 
doping dependence of T^, but can also be observed in ARPES experiments we extended 
previous results for the spectral density of states p(k, u )0 to analyze what happens when k 
crosses the shadow of the FS. In Fig. 4 we plot p(k, u) for k close to the shadow of the FS in 
the neighborhood of the (vr, 0) point. For the overdoped system one sees that the maximum 
II in S(^(k, C(j) of Fig. 3 leads to a corresponding dip in p(k, cj) at ~ +60 meV, where 
spectral weight is suppressed with respect to the normal state. This dip-structure agrees 
well with the experimental finding of an anomalous dip structure at the shadow of the FS 
which is mostly pronounced near the (vr, 0) point with energy Uexp ~ 70 meVSlll. In the 
normal state of the underdoped system the small satellites below the Fermi level, which are 
separated from the main band above the Fermi energy, are the shadow states. They become 
visible below x = 0.13 and increase with decreasing doping. Below Tc the dip structure as 
discussed for x = 0.16 shifts to smaller energies leading to a transfer of spectral weight below 
the Fermi level0. Note that a first experimental indication for this anomalous transfer of 
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spectral weight was recently obtained by Dessau et al.u. 

Summarizing, we presented new results for the superconducting properties and the phase 
diagram of the 2D Hubbard model (Fig. 5). In particular, we obtained an optimal doping as 
a result of the constructive interference between magnetic correlations and ci-wave pairing 
state. Furthermore, we demonstrated for the first time that shadow states exists below 
which are related to dip structures and to an anomalous transfer of spectral weight at the 
shadow of the FS in the spectral density of states. These were observed in photoemission 
experiments indicating the importance of spin fluctuations for the formation of supercon- 
ductivity in the cuprates. 
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particle hole transformation and to exchange momenta k and k + Q. An even better 
agreement with the photoemission experiment might be obtained by using the observed 
topology of the FS in the Bi2212 system. 
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FIGURES 

FIG. 1. Doping dependence of the superconducting state: (a) Tc obtained from A(T = Tc) = 0. 

(b) Enhancement of the effective interaction max(V^(k, 0))/max(V^(k, 0)), 2Ao and inverse quasi 
particle hfetime t~^. (c) Ratio 2Ao/A;bTc. 

FIG. 2. Superconducting gap function A(k, a;) for different frequencies and along different 
paths in the Brillouin zone (see inset). Note that A(k, a;) is normalized by A(k = (tt, 0),a;) and 
that for comparison we plot the simple d-wave gap function A^ ~ cos{kx) — cos{ky). 

FIG. 3. Im S(^(k,a;) at the Fermi surface and at its shadow. Note, the results T = Tc refer to 
the normal state. 

FIG. 4. Density of states p(k, lo) along a path in the Brillouin zone that crosses the shadow of 
the Fermi surface above and below Tc. 

FIG. 5. Schematic phase diagram of the 2D Hubbard model, where the antiferromagnetic 
phase is only plotted for illustration although we get no magnetic phase transition in two dimension. 
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